NASA 

Technical 

Paper 

1948 

December 1981 



NASA 
TP 
1948 
c- 1 



Optimization and PerformaniB || 
Calculation of "Sli 

Dual-Rotation Propellers 

Robert E. Davidson 


' ,'t.O AH COPVC TO r 

,'ArWt,, TECHNICAL LIBRARY, 
.kiirrcANp ,,Ara 1^-. M , , 


rWNSA I 

V 



NASA 

Technical 

Paper 

1948 

1981 


TECH LIBRARY KAFB, NM 




□ □ 





Optimization and Perfomiance 
Calculation of 
Dual-Rotation Propellers 


Robert E. Davidson 

Langley Research Center 
Hampton, Virginia 


lUASA 

National Aeronautics 
and Space Administration 


Scientific and Technical 
Information Branch 


SUMMARY 


An analysis is given which enables the design of dual-rotation propellers. It 
relies on the use of a new tip loss factor deduced from T. Theodorsen's measurements 
coupled with the general methodology of C . N. H. Lock. The analysis eliminates the 
possibility of obtaining an infinite chord as would be found by using the tip loss 
factor advocated by Lock. In addition, it includes the effect of drag in optimizing 
and does not require the averaging of various quantities across the radius in carrying 
out off-design calculations. Thus, a combination of the Lock and Theodorsen formula- 
tions is described and the possibilities are explored. Some values for the new tip 
loss factor are calculated for one advance ratio. The calculation is simple and 
straightforward . 


INTRODUCTION 

The current interest in fuel-efficient air transportation has given rise to a 
number of studies aimed at defining the capabilities of large propeller-driven air- 
craft employing advanced, aerodynamic and engine/propeller concepts. An opportunity 
for designing more efficient wings and propellers is provided by new design tools 
which utilize nonlinear transonic-flow codes and improved materials and structural 
concepts. The development of techniques for achieving more efficient wings has 
received much attention, particularly in the NASA Aircraft Energy Efficiency (ACEE) 
Energy Efficient Transport (EET) Program carried out over the past 4 or 5 years. 
Propeller theory, on the other hand, has been pretty well ignored since the early 
fifties. Much more research needs to be invested in propeller aerodynamics to bring 
propeller design up to the level of sophistication achieved by current wing-design 
methodology . 

There are many papers treating both single- and dual-rotation propellers. 
Analyses pertinent to the present investigation that treat single-rotation propellers 
can be found in references 1 to 6, and in the references therein. Dual-propeller 
analyses can be found in references 1, 3, 7, and 8. In references 1 and 7, the prob- 
lem of minimum induced energy loss is treated at length. Provided in reference 8 is 
an important treatment of dual- rotation-propeller aerodynamics. A treatment of the 
calculus of variations that is well suited to establishing Betz' condition (eq. (1.3) 
in ref. 9) for the optimization of both single- and dual-rotation propellers is given 
in reference 10. 

Overall, it would seem that the theory of single-rotation propellers is slightly 
more advanced than that for dual-rotation propellers. The theories of Lock and 
Theodorsen, when applied to dual-rotation propellers, have their shortcomings. For a 
dual- rotation propeller, any attempt to derive an optimization formula (Betz* con- 
dition) based on the Lock formulation of dual-rotation aerodynamics will be found, in 
practice, to yield a planform which develops infinite chords. The reason is that the 
tip loss factor used by Lock is only known for single-rotation propellers and is 
inadequate for dual-rotation propellers. Theodorsen's analysis does not permit the 
optimization of a dual-rotation propeller with drag considered and, in making off- 
design calculations, requires the use of the mass coefficient which is averaged over 



the radius. The purpose of this paper is to combine the best features of those two 
methods into a single, modified theory which eliminates the shortcomings just 
described. 

A way has been found to determine Lock's tip loss factor for dual-rotation pro- 
pellers from Theodorsen's measurements. In order to accomplish this, it was neces- 
sary to show that the Theodorsen formulation applies to the setup in which one pro- 
peller is behind the other. This new Lock/Theodorsen tip loss factor eliminates the 
infinite chords and permits an analysis which includes drag and eliminates the need 
for an averaged quantity. Thus, a combination of Lock's and Theodorsen's formula- 
tions is described and the possibilities are explored. Some values for the new tip 
loss factor are calculated for one advance ratio. The calculation is simple and 
straightforward . 

SYMBOLS 

A abbreviation variable (eqs. (3)) 

b equation (B23) 

B number of blades in propeller 

c chord of airfoil 

Cj^ drag coefficient of blade 

lift coefficient of blade 

o cr 

Cp power coefficient, fiQ/pn d 

d diameter of propeller 

do drag of blade element 

J advance ratio, V/nd 

k number with same value at all radii to blade sections 

K(x) circulation function used by Theodorsen (refs. 1 and 7) 

K' denotes something to be held stationary, in sense of calculus of variations 

£ reciprocal of lift-drag ratio 

dL lift of blade element 

n revolutions per second 

dP power loss of blade elements 
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p',q',r',s' elementary functions of propeller parameters and functions 
(see eq. (16) ) 

dQ torque on blade elements 

r radial coordinate 

R tip radius 

s solidity of either component, Bc/2TTr 

t-j^ temporary variable (see eq. (7)) 

u,v components of interference velocity of front propeller on back propeller, 

or vice versa (fig. 1) 

V forward speed, or advance velocity (fig. 1) 

w trailing helix displacement velocity 

w = w/V 

interference velocity of either airscrew on itself 
W resultant velocity at blade element (fig. 1) 

Wq W for light loading limit (fig. 1) 

X = r/R 

a angle of attack of two-dimensional airfoil 

3 total induced angle (see eq. (B9) and fig. 1) 

Y self-induced angle (see eqs. (B7) and (B8) and fig. 1) 

^Qy equation (B18) 

n efficiency 

0 blade angle (no load) 

p mass density of air 

a product of solidity and lift coefficient, sCl 

(})q see equations (B9) and (BIO) and figure 1 

c|)gQ resolving scalar (eq. (9) ) 
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Xq tip loss factor (see eq. (Bl) ) , X(^q) 

angular velocity 
Subscripts : 

B back propeller 

C mean value for dual- rotation pair 

F front propeller 

S single propeller 

y,z either front airscrew and back airscrew, respectively, or vice versa in 

equation (B17) 

1 induced 

2 drag 

OPTIMIZATION FORMULA 

Betz' condition concerns the effect of making small changes in the chord or the 
blade angle of a propeller at various radii. It gives a mathematical statement that 
must hold at each radius when the changes have been made in such a way that the pro- 
peller efficiency is highest. For single-rotation propellers, Betz' condition can be 
derived intuitively as was equation (1.3) in reference 9; however, the more analytical 
approach of appendix A is preferable, for dual-rotation propellers, especially if the 
propellers are not assumed to be alike front and back. In this paper, only changes in 
the chord are considered. The lift coefficient is presumed fixed by airfoil 
considerations . 

The equation by which a dual-rotation propeller can be optimized will now be 
developed from Betz' condition for dual rotation. This condition is 


k 


d f^^^F ^ 

da \dr ^ 




( 1 ) 


which is equation (A3) in appendix A. 

Equation (1) has to be true at any radius r along the blade. In particular, 
the constant k is to be the same all along the blade. The exact choice of k 
involves matching the power absorbed by the propeller to the power output of the 
engine. 

The various items in equation (1) will now be selected from appendix B, which 
gives the appropriate part of Lock's theory. Familiarity with this appendix and with 


4 


figure 1 (taken from ref. 8) is assumed. In equation (1) , dPp and dPg represent 
both induced dP-|^ and airfoil-drag dP 2 power losses. The induced power loss for 
both propellers is determined by equation (B33) : 



2AbO^(l + Xq cos 2cj)g) 


( 2 ) 


where 


A = TTPrr^^^^ sec^ (pQ 


O = sC^ 

J-i 

1 /^ = 4Xg Sin c()o 




Y = ba (see eq. (B23) ) 


(3) 


Lock's tip loss factor Xg is of primary importance and more will be said about it 
toward the end of this section. The airfoil-drag power loss for the combination can 
be similarly written, from equation (B25) , 



2AG£ 


(4) 


where 


£ = Cg/CL (5) 

Adding equations (2) and (4) gives the term to be differentiated in the numerator of 
equation (1); thus, 


^i’B\ 

i,dr ^ dr / 


2Aa(t-La + £) 


( 6 ) 


where 


t]_ = b (1 + Xg cos 2 c|)g) 


(7) 


The power input to either propeller is given by equation (B26) as 


Q 



A 


q 

sec c|)g 



( 8 ) 
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where 


4>qQ = sin (pQ + Z cos (J)q 


(9) 


Equations (6) and (8) are substituted into equation (1) and are differentiated with 
respect to a. The result is 


2(2t^a + Z) 

(f)qo/sec (pQ 


( 10 ) 


Equation (10) can be solved for a to obtain 
1 


a = 


2 ’"'i’qO ^0 ■ ^ 

2t, 


( 11 ) 


Substitution of equation (7) into equation (11) yields the optimization formula for 
the dual-rotation propeller. 


a = 


k(f)qo cos (|)q 


- Z 


4 sin (P^[Xo 


+ COS'" 


4*0 


( 12 ) 


in which the third of equations (3) was used. 

Equation (12) is used to determine the blade-chord distribution once the 
distribution over r is known. As noted for equation (1) , equation (12) has to be 
true for all values of r < R, and in particular, k is the same for all r with a 
value that has to be found by trial and error to make the power absorbed by the pro- 
peller equal to the engine power output. The power absorbed is found by graphically 
or numerically integrating equation (8) over the radius for both propellers (i.e., 
the power loss for both propellers is found by multiplying equation (8) by 2) . It is 
helpful to note that the power absorbed should increase with k. Further note that 
instead of regarding lift coefficient as given, the chord distribution could have been 
prescribed and then the optimum lift coefficients determined; the airfoils might then 
be optimized for these lift coefficients. Consider now the function Xq- 


LOCK'S TIP LOSS FACTOR Xq 

At a given advance ratio J, the function Xq is a function of x only. It is 
the ratio of the induced angle of attack with an infinite number of blades to the 
induced angle of attack for whatever number of blades happens to be used. Stated 
another way, it is the average rate of fall of potential, taken around the circle of 
the blade element, divided by the normal derivative of the potential at the vortex 
sheet. In view of these physical meanings, it seems obvious that Xq must be a 
significant link between single- and dual-rotation aerodynamics. 
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It might be thought that the Xq dual-rotation propellers ought to be given 

a new symbol since the significance seems so different from that in single- rotation 
propellers. Strangely, however, the new Xq is still a single-rotation function 
because it is merely the Xq associated with a single-rotation propeller. This 
single-rotation propeller, however, does not have the optimum single-rotation pro- 
peller loading; it has instead the optimum dual-rotation loading which could be 
obtained from Theodorsen's circulation function K(x). This radical change in load- 
ing makes a considerable difference in the function Xq- 


Applicability of Theodorsen's Theory to Dual-Rotation Propellers 

As pointed out in the Introduction of reference 1, Theodorsen's work was based 
on the hypothetical situation in which the dual-rotating propellers are operating in 
the same plane. It was noted that the applicability of Theodorsen's theory to other 
situations requires further confirmation. This matter will now be considered because 
it is pertinent to the determination of Xq dual-rotation propellers operating 

one behind the other on a single axis. 

Observe that, when the propellers are in the hypothetical situation of being in 
the same plane, equations (B32) would be symmetric; thus, 

(^)p = ^ ^0 ‘J’o ~ ^o) 

^ ^0 ^o) 

But, for the combination obtained by adding these two equations, the result would 
again be equation (B33) in appendix B, with no change in equation (2) . Therefore, 
the previous derivation for the optimizing equation would proceed with no change; 
equation (12) for O would be the same whether the two propellers were in the same 
plane or not. So the optimum planform defined by Theodorsen's theory applies without 
reservation when the two propellers are mounted one behind the other. 


Determination of Xq for Dual-Rotation Propellers 

The only quantity in equation (12) that cannot be considered known is Xq^ which 
first appears in equation (2) . The single-rotation values for Xq cannot be used, 
as already noted, but Xq for dual-rotation propellers can be determined from a com- 
parison of the Lock and Theodorsen theories. 

If the physics and mathematics leading to equation (12) are correct, this equa- 
tion must produce, with drag neglected, the same results as are obtained from 
Theodorsen. On page 87 of reference 1 (the equation just before eq. (12)), Crigler 
gives the following equation for light loadings: 




(Mr 


V 


TTndx 


wV K(x) 
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or 


a = w -^(sin (bn) K(x) = s'w 

TTX U 


(13) 


On the other hand, equation (12) with drag neglected is 


a = 


- k(j)qo cos (pQ 


4 sin 


+ cos^ - Sin2 


(14) 


Equating equations (13) and (14) gives 


where 


, Xq 


ws 


(15) 


P' = I V 


q' = 


4 sin (f) 


0 


(16) 


9 p 

r' = cos (})g - sin (})q 


s ' = — K (x) sin 
TTX 


J 


In the last of equations (16), K(x) comes from the electrical measurements of 
Theodorsen (fig. 2). Equation (15) can be solved for Xg : 


X = q s 

® (k/w)p' - q's'r' 


(17) 


An independent derivation has been made of the important equation (17) for 
determining Xg from Theodorsen' s electrical measurements. This derivation deals 
directly with the induced velocities and it is interesting that Betz' condition does 
not play any part, although it did in the earlier derivation. For details, refer to 
appendix C . 
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Everything can be regarded as known in equation (17) except the ratio k/w, 
which will now be found from somewhat tangential considerations. First note that, 
with drag neglected, the efficiency decrement can be determined by using equations (2) 
and (8) , 


2b (1 + Xq cos 2(j)Q)a 2t^a 

2'1’qO ‘t’o ^^qO 


If equation (11) , with £ = 0, is substituted into equation (18) for O, there 
results 

1 - n = ^ k (19) 


dP 

dr 


2f2 


dr 


In addition to equation (19) , another relation can now be found between n 
and w. In figure 3 (which is taken from ref. 1) , q is seen to be practically 
linear with w for light loadings. The slope dq/dw is taken to be 


dq 

dw 


-0.5 


( 20 ) 


which means that 


1 - q 



dw 


0 . 5w 


( 21 ) 


Therefore, from equations (19) and (21) , it follows that 


or 


0.5w 





I = 2.0 

w 




( 22 ) 


With this result, everything is known in equation (17) . 

The tip loss factor Xg defined by Lock can now be calculated from the electri- 
cal measurements of Theodorsen (refs. 1 and 7) by using equation (17) . Some numerical 
values have been calculated and are discussed in "Results" and in appendix D. 
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Limiting Forms of Xq 

The behavior of Xq near x = 0 and 1 will now be investigated. 

X = 0 .- It can be seen that p' = 0 and r' - -1, since (j)Q = Tr/2; and q' = 1/2, 
while s' a- oo (see eq. (16)). Therefore, equation (17) shows Xq = 1 at x = 0. 

This is as it should be, since Xq = 1 would be true in the vicinity of a vortex on 
the axis; this vortex is a feature of the optimum dual-rotation propeller. 

X = 1 .- The functions p' , q' , and r' are all finite (see eq. (16)); but 
s' = 0 at X = 1 because K(x) is zero there. Equation (17) now shows clearly 
that Xq = 0 at x = 1. This is as it should be from the definition of Xq - 

All the equations necessary for the optimization of a dual-rotation propeller, 
with drag considered, have now been obtained. But for complete definition of the 
propeller (i.e., to specify the blade-angle distributions, front and back, in addition 
to the chord distribution) , the equations for the induced angle of attack 3y are 
needed. These equations are picked out of appendix B and given explicitly in the 
section "Dual-Rotation Propeller Calculations," where they are also needed for other 
purposes, like performance calculations. 


Recapitulation 

In the optimization equation (eq. (12)), everything on the right may be regarded 
as something given, except Xq, which is found from equation (17) . In using equa- 
tion (17), equation (22) has to be used. It was encouraging to find that Theodorsen's 
work could be applicable to the setup where one propeller is behind the other. It 
would be possible to determine Xq for all possible dual-rotation propellers, before- 
hand, using equation (17) ; then, in this sense, everything on the right in equa- 
tion (12) would be known. 


DUAL-ROTATION PROPELLER CALCULATIONS 
Optimum Propeller 

In order to calculate the optimum chord distribution, it is only necessary to be 
given the values of J and Cp and the "airfoil data." The "airfoil data" are 
actually preselected in the sense that the airfoil and the value of Cp, Cp, or a 
have been preselected to produce a desirable condition, like ('^l^'^D^ max' each 
radial station along the wing. Then, a can be calculated from equation (12) using 
Xq obtained by the procedure given in appendix D. 

If the airfoil data are limited and therefore erratic, this will be reflected in 
the shape of the blade planform. This could cause the blade to be unacceptably 
erratic, requiring smoothing and implying a questionable smoothing of airfoil data. 
The sensitivity of the airfoil sections to Mach number and Reynolds number may 
restrict the operating range of altitude and flight velocity at which the propeller 
will be optimum. For instance, an operating condition could be envisioned in which 
the airfoils are supposed to be supercritical over much or all of the blade so that 
the planform might take a very special shape. 


10 



The no-load blade angle is given by 


9y. + (Torsional deflection) = <f)Q + 3y + ot 


( 23 ) 


in which the total induced angle is given by equation (B27) (taking account of 

eq. (B31) ) so that 


^y = ^(1 + ^Oy) 


(24) 


and further, from equation (B23) , 


= bad + Coy) 


(25) 


thus , 


The Cqy defined by the equations immediately following equation (B28) ; 


^OF '^0 


Cqf ^0 


^cos^ c{>o - 2 sin^ c{)q) 


(26) 


Then equation (25) can be written out, for use in equation (23) , as 

3p = bad + Cop) 

\ (27) 

3b = bad + Cob) 

Note that equations (27) do not depend on any averaged quantities like the "mass 
coefficient" (compare with p. 86 in ref. 1, relative to interference velocities) . 

The detailed procedure for optimization is given in appendix E. 


Nonoptimum Propeller, Given Propeller, Off-Design Conditions 

Mathematically, these problems involve replacing equation (12) by equation (23) . 
The resulting system, comprising equation (23) , two new induced angle-of-attack 
formulas in place of equations (27) , and the two-dimensional airfoil data (which may 
come from either wind-tunnel test or from airfoil theory) are to be solved by itera- 
tion,' perhaps starting with the assumption that the 3y are zero. Then, an initial 
ay can be calculated from equation (23) and a starting is taken from airfoil 

data. Next, an initial 3y can be found which yields new values of and 

and establishes an iteration loop. 
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A new formula for 3y is needed because the are no longer the same front 

and back (although s is) . From equation (B27) it is seen that the formulas (27) 
for By now become 


or 


By - bs (Cj^y + 


- bs (C]^p + ?OF^LB^ 




Bg bs (Cgg + CggCgg) 


(28) 


A new form of equation (23) should be added to these equations because now ay are 
not the same front and back: 


0y + (Torsional deflection) = c()q + By + 


(29) 


This equation is usually very sensitive because 9y and (J)q are often nearly the 
same. In particular, the torsional deflection may be considerable. 

The nonoptimum propeller calculation is not simple because it is iterative and 
requires the storage of airfoil data. There is also a more fundamental difficulty: 
the circulation function has only been determined for the Theodorsen optimum pro- 
peller. The more the propeller and/or operating condition departs from Theodorsen' s 
optimum, the more questionable this calculation is. However, experience with single- 
rotation propellers indicates that the results of this calculation will hold sur- 
prisingly well. 

A more rigorous, but clearly more difficult, method would include arbitrary 
propeller theory in the loop (as given in refs. 2 and 4). Then, each time a circula- 
tion distribution is obtained, a rigorous induced velocity and By would be found. 
This simpler procedure should be of great value in getting started and often may be 
sufficient without the introduction of arbitrary propeller theory. 

Lock's methodology (ref. 8) appears to make calculations for dual-rotation pro- 
pellers essentially like those for single-rotation propellers. 

See appendix F for detailed procedure of nonoptimum propeller calculations. 


RESULTS 

The function Xq has been calculated for dual-rotation propellers, with 
J = 5.1693 at several values of x, by using the procedure outlined in appendix D. 
There were four blades front and four back. 
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The calculated values of Xq are shown in figure 4 plotted against x. Also 
shown is the corresponding curve for single-rotation propellers. The curves are seen 
to differ considerably. Note that since sin (|)q is a function only of x at a 
given J, Xq can also be plotted against sin (})q as in figure 5. 

Shown in figure 5 (fig. 5 of ref. 6) are the conventional contours of Xq for 
single- rotation propellers. The arrows show how the base points (single rotation) 
are shifted for dual rotation. The arrow points and the base points are calculated 
for the same value of J; hence, the shift is vertical. Observe that by performing 
the calculation of Xq at a sufficient number of values of J, a new set of curves 
can be mapped like the ones for single-rotation propellers. 


DISCUSSION 

The paper is now largely complete. Some isolated topics will now be taken up in 
the light of what has been said in previous sections. 


Question of Infinite Chords 

If equations (16) and (17) are substituted in equation (12) with £ = 0, there 
results 


a 


k d • J. 
sin (pr, K(x) 

k/w ° 


(30) 


but, since 


a - sCj^ 


and 

_ B (c/d) 
7TX 


it follows that 


c^ 

d 




k/w 


J 


K(x) 


sin (J)q 


(31) 


Since all quantities on the right of equation (31) are finite, it is clear that c/d 
is finite. There can be no infinite chords when Xq is determined in the way given 
in this paper. 

Competition Between Single and Dual Rotation 

The equations given herein degenerate easily into single rotation without change 
in form. It is only necessary to see that any term involving ^Qy is to removed. 


13 


The optimization of single- or dual-rotation propellers is a simple calculation; 
therefore, the safest and best method of evaluating single- and dual-rotation pro- 
pellers would seem to be a simple comparison of various complete propeller optimiza- 
tions as to efficiency, weight, and cost, rather than an attempt to discern trends in 
the equations. Cost might be set proportional to the total number of blades. Weight 
might be strongly influenced by J, but it is not too clear because the propellers 
turn slower with increased J and centrifugal stresses are relieved. 

In particular, the two curves in figure 4 labeled "single" and "dual" tell 
nothing of the relative merit of single and dual rotation, because they are both for 
the same advance ratio and the total number of blades for the single-rotation pro- 
peller is only half that for the dual-rotation propeller. The single-rotation pro- 
peller is likely to have twice the number of blades as either of the dual-rotation 
pair and the values of J might differ considerably. 


Comparison of Xq for Single and Dual Rotation 

The most characteristic difference between the two kinds of Xq is that the 
values for single rotation greatly exceed unity for inboard radii while the values for 
dual rotation stay below unity or exceed it only slightly. The reason for this dif- 
ference is that Xq = 1 marks the radius where, for a high-J propeller, the benefit 
from dual rotation completely cancels self-induced losses. 

This cancellation can be seen in the equations for the induced power loss and the 
induced angle of attack. The axial losses do not participate in the dual-rotation 
action and tend to make the cancellation less clear so that, for simplicity, the 
propeller of very high J will be considered. Now, tt/2 and cos 2(J)q -> -i; 

therefore, equation (2) for the induced power loss becomes 



= A2ba2(l - Xq) 


which shows that the induced power loss becomes negative when Xq exceeds unity. 
Clearly, the optimized chords will become large where Xq = 1- Any optimization 
process must have built-in controls for preventing the catastrophe of infinite chords 
at Xq = 1. Theodorsen's electrical-analogy approach bridged all this mathematical 
difficulty and went directly to the ultimate solution. 

These matters can be seen again in equations (26) and (27) for the induced angle 
of attack. There it is seen that the place on the blade where Xq = 1 is where the 
sum of the induced angles of attack of front and back propellers is zero. Again, 
these remarks apply to the very-high-J propeller where the beclouding effect of the 
axial losses is absent. 


Compromised Optimum Propeller 

It seems inevitable that the optimum dual-rotation propeller will be compromised 
because of the awkwardly large inboard chords. In other words, the chords inboard 
may be arbitrarily reduced for a practical reason, like a prohibitively great length 
of propeller shaft needed to accommodate the large inboard chords. 
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Furthermore, it is true that the efficiency is not very sensitive to variations 
of the planform from the optimum, as in wings where the straight-tapered planform is 
almost as good as the elliptical planform. So, why is the optimum given so much 
attention in the literature? Perhaps the answer is that the optimum serves as a ref- 
erence by which compromises can be kept under control. Thus, most of the work of 
aerodynamic optimization may be done from the standpoint of the given propeller 
(app. F) with relatively little attention given to the optimum propeller (app. E) . 

The function Xq is indispensable in calculating the performance of the given 
propeller (app. F) , yet it is determined from measurements for the optimum propeller. 
It is somewhat paradoxical to use off-design calculations to optimize a propeller 
when these off-design calculations make use of a Xq determined from optimum pro- 
peller results. The point is single-rotation experience indicates Xq can be 
"stretched" to provide answers that are useful in off-design conditions, although 
this may not extend as well to dual-rotation propellers. 


CONCLUSIONS 

The Theodorsen and Lock treatments of dual-rotation propellers were combined, 
and it is possible to draw the following conclusions: 

1. The function Xq, tip loss factor, used in the Lock treatment can be deter- 
mined for dual-rotation propellers from Theodorsen' s electrical analogy. Formerly, 
these functions only existed for single rotation and were inadequate for dual 
rotation . 

2. The effect of airfoil drag can be included in the optimization of dual- 
rotation propellers. 

3. Combination of the Lock and Theodorsen treatments enables the off-design per- 
formance of dual-rotation propellers to be estimated without reliance on an averaged 
quantity, such as the mass coefficient advanced by Theodorsen. The mass coefficient 
has only one value for the whole propeller disc. 

4. Conclusion 3 also applies to the calculation of the blade angles of optimum 
dual-rotation propellers. 

5. From the Lock treatment of dual-rotation propellers, it can be shown that the 
optimum planform is the same whether the two propellers are in the same plane or one 
behind the other. 

6. The combination of the Lock and Theodorsen theories appears to enhance both 
of those works . 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
November 23 , 1981 
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APPENDIX A 


CALCULUS OF VARIATIONS APPROACH TO MAXIMUM EFFICIENCY 
The quantity that is to be minimized is the total power loss for both propellers 




+ 


) 


dr 


The quantities that are held constant in the process are the power absorbed for both 
propellers , 



dr 


and 







dr 


The latter two are not summed. They are held constant individually, but it is not 
stated at this point what these constants are. The chords, and hence 0, are not yet 
assumed to be the same front and back. 

Chapter 6 of reference 10 will be followed with particular emphasis on 
sections 6.2 and 6.5. In these sections, there is unfortunate, but probably neces- 
sary, rotation of the meanings of symbols. The independent variable 0 becomes y 
and the dependent variable r becomes x, in paragraph 6.5. But in paragraph 6.2, 

X and y represent dependent variables like 0, and t represents the independent 
variable. 

For K' (in ref. 10 K is not primed) , 


K' 


dPp ^ dPe \ 
dr ^ dr / 


+ ki; 




<3Qp 




dr 


(Al) 


in which the dP and dQ depend on 0. 

Then for Euler's equation (eq- (6-15)) in the reference, 


9k’ d_ ^ 
da ~ dt 30 ^ ° 

There is no in this problem, so the second term is zero. The Euler equation 

becomes 


3k' 

30 


0 
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Since there are two O' s 


{Op and Og) , there are two Euler equations. Thus, 



0 


and 



= 0 


Now the sum of the power losses depends on both Op and Dg, but the power absorbed 
by the front propeller does not depend on Og and vice versa. The two Euler equa- 
tions then become 


) (A2) 


In these equations, the dP and dQ are taken from appendix B, then equa- 
tions (A2) are two relations defining Dp and Og as functions of r. When O is 
the same front and back as in the text, the two equations become one (Qp and Qg are 
the same when Op = Og in the approximation accepted) . Then equations (A2) become 


dPp dP- 


and 


BOpVdr dr 


/dP^r dP^ 


dr j ^ SOpVdr / 


3ap.\dr dr / B 9a 


+ k-D ^ 


3 

3ap,\dr / 


+ 0 


/^p ^ ^g\ ^ ^ ^ 

da \dr dr / dO \dr / 


(A3) 


which, except for the nonessential sign of k, is the same as equation (1) in the text 
and is the desired Betz condition. 
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INTERFERENCE VELOCITY FOR A CLOSE PAIR 
OF CONTRA-ROTATING AIRSCREWS 
by 

C. N. H. Lock 


The British Crown holds the copyright for the report, R & M No. 2084, which is 
reproduced in this appendix with permission of the Controller of Her Britannic 
Majesty's Stationery Office. 
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Interference Velocity for a Close Pair of Contra-rotating 

Airscrews 

By 

C. N. H. Lock, M.A., F.R.Ae.S., 
of the Aerodynamics Division, N.P.L. 


Reports and Memoranda No. 2084 
22nd July 1941 


Summary . — A method is developed of calculating the performance of a pair of contra-rotating airscrews, closely 
analogous to that described in R. & M. 2035^ for a single airscrew. The assumptions made are considered to be 
theoretically justifiable if the interference velocities are so small that their squares and products may be neglected. 
It is hoped to compare calculations by the present method with experimental results. 

The equations have been applied by an approximate single radius method to give the difference in blade setting 
between the front and back airscrews for equal power input ; a comparison is also made between the efficiencies of 
single- and contra-rotating airscrews. 


1. Introduction . — The present note contains equations for a close contra-rotating pair of air- 
screws based on the same assumptions as those of R. & M. 1674^ and 1849^, together with the 
following special assumptions. These assumptions appear to be justifiable when the interference 
velocities are considered as small quantities of 'the first order of which squares and products 
may be neglected. 

(i) The interference velocities at any blade element may be calculated by considering the 
velocity fields of the two airscrews independently and adding the effects. 

(ii) Either airscrew produces its own interference velocity field which so far as it affects the 
airscrew itself is exactly the same as if the other airscrew were absent and includes the usual 
tip loss correction. 

(iii) Added to this is the velocity field of the other airscrew. Since the two are rotating in 
opposite directions, the effect will be periodic and its time average value may be taken to be 
equal to the average value round a circle having a radius of the blade element. 

(iv) In considering the interference of either airscrew on the other, it is necessary to resolve 
the mean interference velocity into axial and rotational components.- 

The average Value round a circle of the axial component interference velocity varies slowdy 
through the airscrew disc. It is therefore reasonable to assume for the axial component for a 
close contra-rotating pair that the effect of either airscrew {y) on the other {z) is equal to the mean 
axial component in the plane of the airscrew disc of {y).* 

The average value round a circle of the rotational component is zero® at any distance in front 
of the airscrew disc and has a constant value at any distance behind, this value being twice 
the mean effective value for the airscrew blade sections. It is therefore assumed as regards the 
rotational component that the effect of the rear airscrew on the forward airscrew is zero ; the 
effect of the forward airscrew on the rear airscrew is equal to twice the mean value of the rotational 
component in the plane of the disc of the forward airscrew with its direction reversed. 


* Varying degrees of closeness might be allowed for empirically by multiphdng «<, by (1 — n) and b}' (1 + fi), where 
/j is a parameter varying from a small value for a close pair to a value near unity for a distant pair. 

(76)1)4) A 
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- 2. Equations of motion will now be written down on the lines of the above assumptions using 
as far as possible the ordinary notation {see Fig. 1). In order to maintain the greatest possible 
degree of generality the equations will be developed to as late a stage as possible on the basis 
of assumptions (i) and (ii) only. Thus either airscrew is subject to its own interference, velocity 
Wy, which is normal to W (Fig. 1) and is given by the usual equation 

Wy = sCtW/4>tsm <f> ; . . . . . . . . . . . . . . . . (1) 

in addition it is subject to the interference velocity of the other airscrew whose axial and rotational 
components will be denoted by u and v. 

The values of u and v according to assumptions (iii) and (iv) may be obtained as follows. The 
mean value Wy of Wy taken round the circle of the blade element is given by the equation 

Wy — sCiWj4 sin <f> 

= XWy , (2) 

and is in the same direction (normal to W) as Wy* Then according to assumptions (iii) and (iv), 
denoting the front and back airscrews by suffices F and B (Fig. 1, b and c), 

Up = Wyg cos 

= ftgWyg cos lf>g, .. .. .. .. .. .. .. .. (3) 

Ug ^ pU)yp COS 4^ p t .. *. .. .. •* .. .. (4) 

Vp = Q , (5) 

Vg = — 2xpWyp sin <f>p (6) 

In what follows the general notation [u, v) will be retained as long as possible. 

The general equations will first of all be obtained in a form convenient for ultimate reduction 
to a first order theory analogous to that of R. & M. 2035® using the following notation. Write 
for either airscrew 

Wy — W tan y , . . . . . . . . . . . . . . . . . . (7) 

(Fig. 1) which by equation (1) implies also 

sCg = 4x sin <j> tan y . . . . . . . . . . . . . , . . . (8) 

Write also (as in R. & M. 1849®) 

= ^0 + • • • • • • • • • • • • • • • • • • (9) 

where 

V = rQ tan 4>o- • • • • • • • ■ • • • • • • • • • • (10) 

Resolving parallel and perpendicular to the direction of W (Fig. 1) for either airscrew^ 

W = rQ sec <^o cos ^ + u sin <j> — "v cos if> , ■ ■ ■ ■ ■ ■ • • • • (H) 

Wy = W tan y = rQ sec <^o sin ^ — u cos <j> — sin <^ , . . . . . . (12) 

where, if assumptions (iii) and (iv) are made, u and v are given by equations (3-6). 

* Strictly speaking the value of <!> corresponding to «, v will differ from that appropriate to Wy but the difference is of 
the second order in WylW and will be ignored. 

t Varying degrees of closeness might be allowed for empirically by multiplying by (1 — /<) and by (1 n) 
where is a parameter varying from a small value for a close pair to a value near unity for a distant pair. 

J For a single airscrew, ^ = y, and the symbol y is not used. 

§ IF is the projection of the broken line C D E A on A B'.Wyis the projection of the reversed line A E D C on B C. 
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For the thrust and torque acting on a blade element we have the usual equations 
dT = N{dL cos <j> — dD sin <f>), 

{l/r)dQ = N{dL sin ^ + dD cos <f>), 

where 

dL \pcW^Cjdr, 
dD — ^pcW^Cpdr, 

so that 

{dT I dr) = nprsW^ (Cl cos <f> — Cp sin <f>), .. .. .. .. .. .. (13) 

(1/r) (dQjdr) = TiprsW^ (Cl sin <l> + Cp cos ^). . . . . . . . . . . . . (14) 

For the total power loss (power input minus thrust power) we have 

DdQ — VdT = NdL (rQ sin 4> — V cos <j>) + NdD {r£2 cos 4> V sin <f>) . 

By the geometry of Fig. 1 it follows that for the induced loss (defined here as the part of the 
power loss depending on the lift of the blade elements), 

dPi = NdL (rQ sin F cos <j>), 

(dPi/dr) = nprsWh'D sec C/, sin /?, . . . . ' . . . . . . . . (15) 

and for the drag loss 

dPi = NdD (rQ cos ^ + F sin <f>), 

[dPijdr) = nprsW^Q sec <f>g Cp cos /S. . . . . . . . . . . . . (16) 

Equations (13-16) are all identical in form with those for a single airscrew. 

Equations (10-16) with (3-6) will be developed into forms analogous to those of R. & M. 1849® 
and R. & M. 1674^ in §7 and §8 respectively. The most practical and useful form is obtained by 
considering and y as small quantities and neglecting squares and products of /? and y for both 
airscrews. The resulting equations analogous to those of R. & M. 2035® are developed in §§3-5. 


3. First Order Theory . — Consider /3, y as small quantities of the first order and wrire 

Vy = 

and 

PQy sin ^Oy ^Oy cos <f>Qy ^Oy ) 

POy cos 4>0y + ’’Oy sin <l>Qy = Coy , 

where either y = F, z = B or y.= B, z = F. 

Thus equation (11) gives for either airscrew, 

Wy = rQy sec + SoyWy, + 0(y®) , 

and substitution in equation (12) gives 

{/S — y)yrOy sec <f>„y = CoyrO, sec <i>o,y, + 0(y=) , 

for either y = F, z = B or y = B, z — F. On the basis of equations (3-6) we have, 
~ ^OB COS (f>ag 0{y) , 

*’oF = 0(y) , 

PoB ~ ^OF COS <j>op -j“ ^(y) I 
VoB = — sin <f>op + 0(y) . 

( 76104 ) 




(17) 

(18) 

(19) 

( 20 ) 

( 21 ) 
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Since 

rQy tan < 1 ^ 05 ,=^^ = ^Q, tan <j>o, , 
equation ( 20 ) may be written in the form 
— y)y= ^0y^y.Y. + O(y^) , 
where A,, = (sin ^oy/sin <!>„,). 

Also from ( 8 ) for either airscrew, 
y = bsCi, + 0(y®) , 

where 

l/b = 4»o sin ^ 0 *. 

so that sCl is of the same order as y. 


( 22 ) 


(23) 


If Cl is given for both airscrews, equations (23) determine y and equations (20, 18 and 21) 
determine jS for both airscrews. Then equations (15) and (16) in the form 

[dPijdr) = nprs.CiT^ii^ sec® i^o, /S + 0(y®) , . . . . . . . . . . . . (24) 

[dPijdr) = Tiprs sec® <^o + 0 ( 7 ®) . . . . . . . . . . . . (25) 

give the power losses of either airscrew. In general it is convenient to consider sC^ as a small 
quantity of order y® so that both dPxjdr and dP^jdr are of order 7 ®. To the first order the power 
input to either airscrew is 

Q{dQjdr) = nprsr^Q^ sec® <f>^ (Cl sin + Cd cos ^ 0 ) + 0(y^) . . . . . . . (26) 

The further development analogous to that of R. & M. 2035® required to determine Cl for 
either airscrew for given blade angle setting is given in §5, but it is convenient first to consider 
the application of equations (24-26) to determine explicitly the power input and power wastage 
to the first order, for given Cl, for the particular case of equal rotational speed and power input 
for the two airscrews. 


4. Special Case. Equal Rotational Speed and Power Input. — Equal Rotational Speed. — It 
follows from equations (10) and (23) that equal rotational speed of the two airscrews implies 
equal values of <f>o, }<o and b so that is unity. Equation (22) then gives 

^y = Yy+ UyY. + 0(y®) , (27) 

SO that from (24) 

(dPJdr)y = 7tpr*Q^ sec® <^o (sCi)^ (y^ + Coy y,) + 0(y^), (28) 

and using equations (21) 

foF = cos® <l >0 + 0(y), Cob = (cos® <f>o — 2 sin® 4>„) + 0(y) 

and 

(dP^ldr)p = npr (rQ sec <^o)’ (sCJ^Iy^ + x,,yg cos® + 0(7®), (29) 

{dPJdr)g = npr (rQ sec <j>oY (sQs{yB + ^cos® — 2 sin® <f>„)} + 0(t®) . . . (30) 


Equal Rotational Speed and Power Input. — Equation (26) shows that equal power input to 
the blade element at radius r combined with equal rotational speed implies that 

(sCl)f — {sCs)b = 0 ( 7 ®) 

and . . . . . . . . . . . . . . (31) 

Yf —Yb = 0 ( 7 ") ; 

equations (29) and (30) then become 

(dPJdr)p = npr (rQ sec <f>oY SC 2.7 (1 + r-o cos® <f>o) + 0(y^), 

(dPijdr)B = npr (rQ sec ^oY ^^lY (1 + cos® <f>o — 2xo sin* ^„) + 0 ( 7 ®) (32) 


♦ R. & M. 2035®, equation (10). 
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For the combination of two airscrews 

(dPijdr)c = npr {rQ sec <f>af 2 sClY (1 + cos 2<f>o) + O(y^) . . . . . . . (33) 

Equations (31-33) and (25) transformed into equations for the coefficient p^, pa of induced 
drag power loss, analogous to equations (31) and (33) of R. & M. 2035®, may be used to calculate 
the power loss grading for all radii for a given distribution of sC/, (equal for the two airscrews) ; 
the corresponding blade angle distribution may be obtained from §5. The power input grading 
(torqne grading) may be obtained from equation (26) or more accurately (as in R. & ]\1. 2035^) 
from equation (14) using the more accurate value of W obtained below in §6. In the latter 
case the power input will not be exactly equal for the two airscrews if the values of sC^ are equal. 
The second order difference in sC^, required to make the power inputs equal to the second order 
is determined in §6. Or, the performance for a given blade angle distribution may be deduced 
from the equations of §5 ; the blade angles at standard radius (0-7) might be adjusted to give 
equal power input at that radius. 

Example. — For the purpose of illustration equations (33), (25) and (26) have been used' to 
calculate the partial efficiency for a section at standard radius (0-7) for equaTrotational speed 
and power input. The formulae (deducible from equations (31-33), (25) and (26)) are 

j yCi (1 + Up cos® ^o) -T Cp 

~ ~ cos <j>p {Cl sin ^0 + Cfl cos 4) ’ 


(34) 


1 




yCl (1 + «'o cos® 4>p 
cos (f>p {Cl sin 


— 2xq sin® <f>o) + Cc 
Po Cd cos 0o) 


with 


j jyCL (1 + >^0 cos 2<f>„) + Cjj 

“■ cos ^0 {Cl sin <l>p + Cd cos <f>p) ’ 


(35) 


(36) 


Y = sClI{4xo sin (f>o) = BsCl (37) 

In Fig. 2 values of (1 — rjc) a-re plotted for a range of values of J for (1) a pair of contra-rotating 
two bladers and (2) a pair of contra-rotating three-bladers and for the following values of s. Cl 
and Co : — 


s = 0-090, 


Cl = 0-56, 


Co - 0-017. 


The values of s and Cl are those at radius 0-7 for airscrew R in R. & M. 2021*, while the value 
of Cx) is adjusted to give a partial efficiency for this radius equal to the calculated efficiency 
(0-878) for the whole airscrew. The calculations correspond therefore, to a power input to each 
airscrew of 2,000 h.p. at 450 m.p.h. equal to that assumed in R. & M. 2021* (a total of 4,000 h.p. 
for the two airscrews) for the same diameter, rotational speed and height. They were made 
for a range of values of J from 1-27 to 4-54.* They are compared with the corresponding 
efficiency figures for a single airscrew of double (the same total) number of blades and solidity 
and also with airscrews having the same number of blades as one of the contra-rotating pairs 
and the same total solidity. The equation corresponding to (36) for a single airscrew is 

1 yCl d~ C o ^ 

~ cos (f)p {Cl sin ^0 + Co cos ^„) ‘ ‘ * ’ ’ ' • - 1 / 

with (37) in which it must be remembered that values of Xp and s must be used, appropriate to 
the total number of blades and solidity. Thus the value of s for the single propeller has twice 
the value for the corresponding contra-rotating pair, and so the value of y in (^) would be double 
that in (36) apart from the change in Xp due to doubling the number of blades. 

The results of Fig. 2 show that for the present case the increase of efficiency as between the 
2-bladers (contra-rotating) and the 4-bladers (single-rotating) varies from 1-0 per cent, to 4-6 
per cent., and the increase as between the 3-'bladers (contra-rotating) and the 6-bladers (single- 
rotating) varies from 1 - 7 per cent, to 4-8 per cent, for the particular values of s, Cl and Co chosen. 

* The actual efficiency figures for the highest values of J would in practice be reduced by the increase of Cp due to 
increased compressibiKty effect. 
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5. The Relation between sC^ and Blade Angle d to the First Order, for the General Case.— -This 
may be obtained by a similar method to that of R. & M. 2035®, §3, as follows : — 

Write 

d <f>Q -|- S = O y •• •• •• *• •• •• •• •• .. •• 

and 

asCi = a -4- e 




(40) 


where a and s define the (straight line) lift curve as in R. & M. 2035®, equation (11), and are, 
in general, functions of the Mach number. 


Comparison of (40) and (23) gives 

60 = 6/S + ay + 0{y^) , . . - . . (41) 


which with (22) determines 0p, as functions of yp, y^ and so of (sCi)p and (sCi)g in the form 
0, = yy + + O(y^) (42) 


withy =F, z = B or y = B, z — F. Using the relation = I, the pair of equations 

represented by (42) may then be solved for y^, y, in the form 

(sCJj, = (y/&)y = {(a -f- 6), 0^ — b,Coy^y, 6>,}/{(« + b)p [a + b)g — bpbgtop^oB) + 0{y^) , (43) 

which reduces to equation (13) of R. & M. 2035® on putting Cof = fos = 0. In this pair of 
equations, using (18) and (21) we have 

CoF ^OB ^OB CDS ”1“ ^(x) J • • ■ • ’ ■ • * • • • • (44) 

= «0F (cos <f>0B COS <l> 0 F — 2 sin sin + 0{y) , (45) 

and 

^FB = 1/^BF = sin 9 SoF/sin (46) 


Special Case. For Equal Rotational Speed, using the results of §4 equation (43) becomes 
(sC£)p = y pjb 

= {{a + b) ©p — x^b cos^ if>„©g}/[{a + by — cos^ 

(cos^ 9^0 — 2 sin^ 4,0)} -f 0(y*) , 

(sCi)^ = y gib 

= {{a + b) ©g — x„b (cos* '^0 — 2 sin* 4>o) ®f)/{(« + bf — x^¥ cos* ^0 

(cos* <Ao — 2 sin* .^ 0 )} + Q(y“) . • • (47) 

For equal rotational speed and equal power input to the blade element at radius r equation (42) 
becomes (using 31) 



and so 

©p-Bg==y {Up - Ub) + 0(y“) 
= 2yxg sin* 4>q + 0{y^) 


= ^sC[_ sin 4>o + 0(y'‘) . . . . . . . . . . . . . . . (48) 

This value is plotted against J in Fig. 3 for the values of sCg used in §4 and varies from 0-7 deg. 
to 1 - 3 deg. over the range of J considered. 
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6. Values of W and Q[dQjdr) to the Second Order in y . — The value of W to the second order can 
be obtained from equation (11) in the form 

= r£2^sec + fo/f2,sec 

= rQy sec {1 + + 0(y*) . . . . . . . . . . . (49) 

The expressions for dTjdr and dQJdr involve the factors sin <t>, cos 4> which may be written, 

sin ^ = sin ^ (1 + /3 cot <f>o) + 0(y==) (50) 

cos <f> = cos ^0 (1 — /5 tan ^o) + 0{y^) (51) 

with 

= Yy+ Uy^ytY, + 0{v“) , (52) 

from (22). Equation (49) then gives 

Wf sin 4 ,^ = r^Df tan <f>^^ sec ^oy {1 + [2foy + foy cot 

+ y, cot 4>.y) + o(y^) , • . • • (53) 

Wf cos 4>^ = r^Qf sec <f>^^ (1 + [2f„, — foy tan <^„^] l^.,y, — y^ tan + 0{y^) ; . . (54) 

or 

Wf sin <!>y = r^Qf sec’' ^ov {sin <^oy + i^y.y, [^oy (3 — cos 2<^,,) — sin 2<^oy] 

+ yj, cos 0oy} + 0{y^) , . . . . (55) 

Wf cos 4>y = r^Qf sec" <^oy (cos <^oy + Ooy sin 2<^oy — I’oy (3 + cos 2<^o>)] 

— Yy sin 4>o,} + 0(y“) . . . . . (56) 

In evaluating Q[dQjdr) and VidTjdr) it is reasonable to consider CojC^, as before, as a small 
quantity of the same order as y, and to write 

Q{dQldr) = npr^QW^ sin <f> sC^ (1 + {CojCf) cot <f>a) + 0{y^) , . . . . . . (57) 

V {dTjdr) = npr^Q tan <j>a cos <j> sCi (1 — {Co/Cf) tan ^ 0 } + <5(y®) . . . . . (58) 

In these expressions IT" sin <f>, IF" cos 4>, a-re given by equations (53-56) in which y is given by 
y = hsC^, 

so that the torque and thrust power loss grading may be evaluated as far as terms of order y" 
if the value of sC^ is known to this order for each airscrew. Equations (43-46) give the values 
of sC[_ for each airscrew in terms of the blade angle settings. 

Strictly speaking, equations (23) and (40) are only correct to the first order in y and a, but 
it was suggested in R. & M. 2035® that in practice the curves of against a in the unstaUed 
range and of sC^ against y over a considerable range of large values of J, are straight lines to a 
higher order of approximation. The additional order of accuracy would then apply to equations 
(43) since they are deducible from (23) and (40) by linear transformations ; the values of sC^ 
deduced from (43) for given blade angles would then be sufficiently accurate when substituted 
in (57) and (58) to give values of thrust and torque power correct as far as terms in y". In any 
case the value of power loss given by taking the difference between power input deduced from 
(57) and useful power deduced from (58), will be consistent with (24) and (25) and correct to the 
same order as the latter equations. 
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Case of Equal Revolutions to Second Orrfer.— Substitution of 

~ Mub ~ ^0 COS <^o , 

^OB — 2?tfo sin (f>Q 
in (53-56), gives 

Wf- sin 4>f = tan <f>o sec 4>o {1 + 7f cot 4>o 

+ XoVb (cot 4>a + sin <l >0 cos <^o)} + 0{y^) , . . (59) 

W/ sin <f>B = tan <f>o sec <f>a {1 + Yb cot <f>o 

+ x^Yp (cot ^0 + 3 sin </>o cos .^o)} + 0{y^) , . . (60) 

Wp^ cos <l>p == sec <^o {1 — Yf tan 4>o + Xgyg sin <j>o cos ^o} + 0{y‘) , . . . . (61) 

Wg’‘ cos <j>B = sec <j>g (1 — yg tan <j>g + «o7f (2 tan + 3 sin <j>g cos <^o)} + 0{y^) . (62) 

Equal Power Input to Second Order . — It is evident that the difference C^p — C^^will be of order 
y^ and it is therefore reasonable to assume that Cgp — Cgg is of order y®. The condition of equal 
power input will therefore be taken as . 

(sC^IT^ sin <i>)p = [sCJV^ sin<i>)g + 0{y^) (63) 

Condition (63) may be satisfied by writing == y in (59-62), since this is true to the first 

order, and putting 

(sCJf = sC^ (1 + «oX sin ^0 cos ^o) , (64) 

(sCj)g = sCg (1 — «oX sin ^0 cos ?io) • • • • (65) 

in (63), where sC^ is a mean value between the two airscrews. The final expressions for the 
thrust and torque grading will be for either airscrew, 

Q{dQjdr) — npr*Q^ tan 4>o sec ^o{sCi{l + y [cot 

4- ^0 (cot pg + 2 sin <f>g cos <^o)]} + sCg cot 4>g} ■, . . (66) 

and for the front and back airscrews separately, 

V{dTjdr)p = Ttpr^Q^ tan <^o sec ^o(sQ{l + y [— tan 4>o 

+ 2xg sin ^0 cos sCg tan <^„) , . . , . (67) 

V{dTldr)g = npr*£2^ tan <j>g sec <l>o(sC^{l + y [— tan <j>g 

+ 2xg (tan <j)g + sin <f>g cos <^o)j} — sCg tan ^o) • (68) 

The difference between these expressions for torque and thrust power agrees with the first order 
value of power loss given in (31-33). Expressions for the blade angle to the second order could 
be deduced from §7, equation (76) below, but would be rather complicated. 

7. Exact Transformation of Equations (11) and (12) into a Form Analogous to the Equations 
of R. & M. 18492.— Write 

Uy =Py Wl, , 

Vy = Vy Wi, , 

where either y = F,z = B, or y — B, 
equations (3-6), of (j>, only). 

Write 

p, sin <j> — V cos (f> = ^ , 
p cos ^ + »> sin ^ f . 


(tjy; 

z = F, and Py, Vy are functions of <j>y, 0, (according to 


• • (70) 
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These definitions are analogous to the first order definitions of (17) and (18). Write also 
rQ sec <j>^ cos B — C , 

. . . . . . . . . . . . . . . . . . (71) 

• rQ sec <f>a sin p — D ; \ 
also by (10) 

{rQ tan ^„)f = tan 4,g)g . 

Then equations (11) and (12) become 

W, = C,+ (72) 

w>iy = Wy tan yy== Dy— CyW^, . (73) 


The pair of equations, 

Wly+CyWl, = Dy 
C.Wyy +W„=D, , 

may be solved giving 

Wry = {Dy - ^yD,)/{l - CvC.) 

Then (72) gives 

Wy = {Cy {I- ^y:,) + ^yD, ~ ^ y^ .D y) / [l ~ f,f.) . 


(74) 

(75) 

(76) 


and substitution in (73), using (71), gives 


tan Yy = tan /3^(l — ^y ^')/ 


1 - f A rh - f.?. 


Dy 

Cy 


= tan 



. sin <f,y sin 

’’ sin <t>, sin /S// 1 '' 


sin <f>t« sin /?, , , 

sin <^ 0 , cos Py ~ 


The two equations (77) determine Yf> Yb so in virtue of (8) (sCJ^, and (sCJ^ as functions of 
<I>P, 4>b only, 4>of< 4‘ob being known. Since equations (3)-(6) are only claimed to be correct to the 
first order, the advantage of the present equations over first order equations is doubtful. 

It would be possible to plot (sCJ^? against <f>p, giving for each J a series of curves for various 
values of <f>g and similarly for (sCg)g, (Fig. 4). It would then be necessary to determine inter- 
sections with (sCJ jr against and (sCjg against curves giving consistent values of and <f>g 
and this could be done by a very rapid successive approximation between the two figures. This 
represents the analogue of the use of Chart I in R. & M. 1849^. 


8. Equations of the type of R. & M. 1674k — From Fig. 1. 
AC = Wr cosec y . 

Resolving parallel to AF, we have 

AC cos = yQ — V , 

giving 

Wr = {rD — v) sin y sec {<j> — y) 

= [rQ — v) tan y sec ^/(l -h tan y tan <f>) , 
with tan y = sCJ4x sin <f> . 


(78) 
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For the front airscrew v = 0 and the equation becomes identical with equation (8) of R. & M. 
1674^. For the back airscrew vjrQ is of order y and might be calculated by writing <x.p = a.g. 

V is then given by (Fig. 1) 

V = rQ tan - GC — CD - HK 

= rQ tan <f> — Wi sec ^ — u — v tan ^ . . . . . . . . . . . (79) 

The most convenient form for W is (Fig. 1) 

W = HA - GB - HG 

= {rQ — v) sec (f> — tan <fi , . . . . . , , . . . . . (80) 

which is identical with equation (2) of R. & M. 1674^ for v --- 0. 

The equations (78-80) may be transformed so as to involve non-dimensional coefficients only, 
by dividing by convenient multiples of RQp, RQp. 

The solution of the equations by the methods of R. & M. 1674^ is straightforward apart from 
the occurrence of the term involving v in equation (78) for the back airscrew. A suitable series 
of values >of the blade incidence a is first chosen for both airscrews for a series of standard values 
of the radius. Values of Cl, Cd for either screw are supposed known as function of a, and </> is 
deduced from the equation 

<f> — 6 — a . 

I 

Equations (78), (79), (80), (14), (15) and (16) then determine in succession values of Wi, W, V, 
Q{dQldr), dPildr, dP^jdr (or of suitable coefficients of them) for both airscrews. In evaluating 
the term v in equation (78) it should be sufficiently accurate to write = a^. It is finally 
necessary to plot values of V or of its coefficients Jp and Jp and of Q(dQldr), dP^jdr, dP^jdr or 
their coefficients against a, so as to deduce values of the thrust and power coefficients for the same 
values of F at all radii before plotting against the radius r* and integrating to obtain the power 
input and power loss on the whole airscrew. 

9. Recapitulation. — §1. Of the four basic assumptions as set out in §1, the first two are con- 
sidered to be of general application to an airscrew, subject to any type of external interference. 
The development of the equations is carried as far as possible without reference to the third 
and fourth assumptions and these may require further empirical modification and would in fact 
be modified as a result of increasing the distance between the two airscrews or varying their 
diameters, etc. 

§2. Equations are given of the most general form consistent with assumptions (i) and (ii) and 
determine the total velocity W and the interference velocity of either screw on itself, in terms of 
rQ, <f>, and (u, v) the components of the interference velocity of the second screw ; also for the 
thrust, torque and power loss grading in terms of W, Cl, Co and <j>. 

§3.' In this section squares and higher powers of the interference velocity ratio are neglected. 
This is probably not a serious limitation since it is very doubtful whether the original assumptions 
hold beyond the first order in the interference velocities. Explicit equations are given for 
(dP[ldr), (dQldr), (dP^ldr) to the first order. 

§4. The equations of §3 are applied to the particular case of equal rotational speed and equal 
power input. Explicit equations are given for the partial efficiency at a given radius and for 
the induced loss for front and back airscrews separately. 

* Coef&cients of the t 3 'pe Ai. Pc 2 plotted against = (r/f?)®. 
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§5. This section gives first order results for given blade angles and also the first order difference 
of blade angle between front and back airscrews for the case of equal angular velocity and power 
input. This completes the formulae necessary to obtain the numerical results given in the 
present note. 

§6. Values of W, {dQjdr) and {dTjdr) are given to the second order for known values of C^. 
Difference of between the two airscrews is determined to the second order for equal revolutions 
and power input. The resulting value of the difference between the thrust power and torque 
power checks with the first order estimation of power loss in §3. 

§7. In this section equations are obtained analogous to those on which the charts of R. & M. 
1849^ are based. 

§8. In this section equations analogous to those of R. & M. 1674^ are developed which could 
be used in the absence of charts to calculate the exact performance of an airscrew on the basis 
of assumption (i)-(iv). 


a 

b 

B suffix 
c 

C suffix 
C, D 

dD 
F sufi&x 

J 

dL 
N 
dP, 
dP^ 
dQ 
r. R 
s 

S suffix 
dT 

U, V 

V 

W 

Wo 


LIST OF SYMBOLS 

Reciprocal of slope of lift curve. Equation (40). 

Equation (23). 

For “ back airscrew ”. 

Blade chord. 

Mean value for contra-rotating pair of airscrews. 

Equation (71). 

Lift and drag coefficients of blade element. 

Drag of blade element. 

For “front airscrew". 
nVjRQ. 

Lift of blade element. 

Number of blades of either component. 

Induced power loss for blade elements. 

Drag power loss for blade elements. 

Torque on blade elements. 

Radius of blade element, tip radius. 

Solidity (= Ncj2jir) of either component. 

For single airscrew. 

Thrust on blade elements. 

Components of interference velocity of front airscrew on back airscrew or 
vice versa (Fig. 1). 

Forward speed (Fig. 1). 

Resultant velocity at blade element (Fig. 1). 

Fig. 1. 
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List of Symbols — continued. 


K'l 

Wi 

y, z suffices 

0 suffices 
a 

P 

V 

e 

C, fo 
n 
e 
0 
X 


Interference velocity of either airscrew on itself (equation (1)). 

Equation (2). 

Denoting either front airscrew and back airscrew respectively, or vice versa 
(equation (17)). 

Indicating limiting value for zero lift. 

Blade incidence. 

Fig. 1 ; equation (9). 

Fig. 1 ; equations (7) and (8). 

Zero lift angle ; equation (39). 

Equations (70), (18). 

Efficiency. 

Blade angle ; equation (39). 

Equation (39). 

Tip loss factor ; equation (1). x^ is written for >c{<f>a)- 
Equation (22). 


>Equations (69), (17). 

V, VoJ 

f, fo Equations (70), (18). 

(}>, <l>o Fig. 1. Equations (9), (10). 

D Angular velocity in radians per second. 


Note. 0{y^). The notation used in equation (19) e/c. The statement 

F{v) = F,{y) + O(y-) 

implies that F{y) can be expanded in powers of y in the form 
F{y) = /o + y/j + y */2 + • • • 

and that 

Fi{y) =/o + yfi + y% + • ■ • + r”"'/.-! • 


1 Lock and Yeatman . . 

2 Lock 

3 Lock, Pankhurst and Conn 
, 4 Pankhurst and Fowler 

5 Lock 
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Fig. 3. Values of (0, — & b ) (independent of number of blades) Calculated for Same Conditions as Fip. 2 (for 
Equal Rotational Speed). The Curve also Represents the Differet\ce of Blade Angles (0, — 0,) provided that 

the Zero Lift Angles are the Same for Both Airscrews. 


Fronb akirscrew 


&9Ck gkir^crew 



Fig. 4. Charts Analogous to Chart I of R. & M. 1849’“. 


(7fiJ6<) Wt. 10/7116 6/47 Hw. G.377/1. 
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INDEPENDENT DERIVATION OF LOCK'S TIP LOSS FACTOR FORMULA 

In the text, Xq was found by making O, or optimum planforms, the same between 
Lock and Theodorsen. Now, however, the induced angles B will be made the same 
between the two formulations . 

The first thing will be to specialize Lock's theory to apply where the pro- 
pellers are in the same plane . (The determination of Xq is the same whether they 
are or are not in the same plane, as was seen in the text.) This specialization 
makes the two equations for ^gy> between equations (28) and (29) in appendix B, 
become symmetric so that they are one 


Co = Xq(cos^ (pQ - sin^ (})q) = r'Xf 


(Cl) 


Then the induced angle of attack, equation (B27) , becomes 


6 = b(l + Cq)o = b(l + r'Xo)cr 


(C2) 


for either propeller, where 

^ — -■ — (C3 

4Xg sin 4>o 

from equation (B23) . 

Now B is the induced angle of attack at the propeller plane, not the far wake. 
The induced angle of attack from Theodorsen, therefore, must be multiplied by 1/2. 
Also, the displacement velocity w multiplied by cos ())q gives the velocity of the 
vortex sheet normal to itself. If this is divided by Wg it is converted to an 
induced angle. Thus, the equating of the induced angle, between Lock and Theodorsen, 
yields 


2 Wq 


cos 



(C4) 


with the right side given by equation (C2) . 

Next, Theodorsen' s circulation function will be introduced, as in the text just 
before equation (13) , 


sCt-W,. = aWn = g— wV K(x) 

L 0 0 TTndx 


(C5) 
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or 

=• = sin d)„ K(x) 

W TTX ^0 

where 

w = w/V 


Now equation (C4) can be written, by using equation (C2) , 

2 ^0 ^ 3 ^ A = ^ ^'^0^ £ 

Wq w Vw ^ w 


By using equation (C3) , equation (C8) becomes 


cos 4) 


0 


4 sin <j) 


(X( 


+ r' 


0 \^0 
V 


— K(x) sin (pQ 


(C6) 


(C7) 


(C8) 


Substituting V/Wq = sin c|)q (fig. 1 of app. B) and multiplying by 4 sin (|)q result 
in 


sin 2<t)Q 




K(x) 


Solving for Xq gives 


^0 


TTX 


JK(x) 


sin 2(})n “ 


(C9) 


(CIO) 


Now, the corresponding formula in the text is equation (17) , which is 


><0 


1 

kp' 
q' ws ' 


r ' 


(Cll) 


Compare equations (CIO) and (Cll), which should be identical. They will be identical 
if 


kp' 
q' ws ' 


TTx 

JK(x) 


sin 


2(pQ 


(C12) 
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Substitute the following from equations (16) in the text, into equation (C12) : 


P’ " i ‘^qO '*^0 "" i '*’0 '*^0 


4 sin (J)q 


s' = — K(x) sin d), 
TTX 


= 0 ) 




(C13) 


Then the right side of equation (C12) is 


TTX 1 

— r 7 — sin d)n sin 2(j)n = -~r sin sin 2(^ 

JK(x) sin c|)q ^0 s' 


0 


(C14) 


and the left side is 


1 

k 2 


— sin (jjQ cos cf>Q 


w 


k 1 . A 1 • A A 

^ — 2 sin cPq - sin (})q cos cJ)q 


4 sin (pQ 


k 1 ^0 

w 2 s' 


;in 2 (|)q 


(CIS) 


Equating equations (C14) and (CIS) gives 


1 ^ = 1 (C16) 

2 ^ 


which must be true in order for equations (CIO) and (Cll) to be the same. But equa- 
tion (C16) is in fact true, as may be seen from equations (22) in the text. 

Hence it is concluded that the two formulas, (CIO) and (Cll), are indeed 
equivalent, and it appears to be verified that the method given in the present paper 
for determining Xg, appendix D, is correct and consistent with Theodorsen and Lock. 
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CALCULATION OF Xq FOR DUAL- ROTATION PROPELLER AT A 
GIVEN VALUE FOR J 

1. Choose a set of values of x corresponding to stations or sections on the 
propellers at various radii. These values of x might well be those found in fig- 
ure 2. The entire calculation is performed for each value of x. 

2. Calculate the advance angle (J)q from 

4>o = 90° - tan“^ x^ 

3. For the values of x chosen in step 1, read a set of K(x) from the circula- 
tion function curves (e.g., fig. 2) for the number of blades being considered (four 
front and four back for the purposes of this paper.) 

4. Calculate p' , q' , r', and s' from equations (16). Note that, since ^ 
is taken to be zero in this calculation, (J)qQ = sin (|>q. 

5. Calculate Xg from equation (17). (The value of k/w was determined to 

be 2.0 in the discussion following equation (17). (See eqs . (22).) 
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OPTIMIZATION OF DUAL-ROTATION PROPELLER WITH DRAG CONSIDERED 

1. A value for J and a value for Cp must be given to start. 

2. Tabulate Xq from step 5 in appendix D. 

3. Tabulate the preselected C^, £, and a for each x. For instance, these 

might be for maximum lift-drag ratio all along the blade. 

4. Calculate (f)qQ from equation (9). 

5. Calculate a from equation (12) . In this step, a trial value for k is 
needed. Equation (19) is an aid in guessing this. 

6. Integrate dCp/dx over the range of x, from body to tip, and compare the 
resulting Cp to the given Cp, which must match. 


dCp 

dx 


tt4 

2 


sec 



(for both propellers) 


This equation can be derived from equation (8). Set = r^fi^ sec^ c})q. Change r 
to xR and use Cp = nQ/(pn^d5). 

7. Calculate 

c_ _ 7T xO 

d “ B Cl 

8. Find 3p and 3g from equations (27) . 

9. Calculate the blade-angle distribution 0 from equation (23). 

9p + (Torsional deflection) = (})q + 3p + oi 
0g + (Torsional deflection) = (J)q + 3g + a 
Items 7 and 9 define the propeller. 
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CALCULATING PERFORMANCE OF GIVEN DUAL- ROTATION PROPELLER 
AT OFF-DESIGN CONDITIONS 


The problem may be stated as: given the propeller at some blade-angle setting 
and the value of J, determine the torque and thrust, or determine the power 
absorption and efficiency. 

The existence of initial distributions of (Cl)p and and hence of Op 

and Ug can be assumed. This may be a guess, or else they can be calculated immedi- 
ately by assuming the induced angle of attack is zero. 

1. Follow the steps in appendix D to obtain the function Xq applicable to the 
given J. 

2. Calculate cfjgQ from equation (9) (the same as step 4 in app . E) . 

3. Tabulate or store airfoil data so that Cg and £ can be found for any value 
of a. (Note that this is not the same kind of airfoil data as in step 3 of app. E, 
which is preselected airfoil data.) 

4. Tabulate 0p, Op, and s from propeller-geometry information and the given 
blade-angle setting. 

5. Calculate Bp and 3g from equations (28). Initially, these might be set 
equal to zero. 

6. Solve equation (29) for and ttg. 

7. Find new (Cp,)p and (CL)g from airfoil data using Ctp and from 

step 6. Also find new Op and Ug. 

After completing step 7, it is possible to return to step 5 and loop through to 
step 7 repeatedly until (Cg)p and (Cg)g converge on final values. At step 7 on 
the last loop, find £p and £g in addition to the lift coefficients. Then step 6 
in appendix E shows how to get the power input. The efficiency is found as in equa- 
tion (18) , but by using equation (6) instead of equation (2) . Note that in place 
of O there are Op and Og. Therefore, the power formula in step 6, appendix E, 
has to be appropriately modified. The same applies to equation (6) in the efficiency 
calculation. This would conclude the calculation if arbitrary propeller theory were 
not to be added. 

If arbitrary theory is to be introduced, the return to step 5/ after completing 
step 7 on the last loop, would call for arbitrary propeller theory rather than equa- 
tions (28) . More specifically, the circulation can be easily found from (Cl)p and 
(Cg)Bf which are found in step 7. From the circulation, the single-rotation part of 
the induced angle of attack Yy can be found, by arbitrary propeller theory, so 
that By now contains arbitrary propeller theory. The tip loss factor Xq is no 
longer used. The loop continues on to step 7. This should produce a rigorous 
lifting-line performance calculation, but the introduction of arbitrary propeller 
theory has to be regarded as a considerable escalation of computational effort. 
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Plane of rotation 
(a) Either airscrew. 
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(b) Front airscrew. 
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(c) Back airscrew. 


Figure 1.- Interference velocity components for pair of 
dual-rotation propellers. 
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Figure 2.- Circulation function K(x) for dual-rotation 
propellers; four blades front and four blades back 
(ref. 1) . 





Figure 4.- Tip loss factors for single- and dual-rotation propellers 
plotted against x; J = 5.1693, four blades single, eight blades 



Figure 5.- Typical presentation of Xg for single rotation for 
four blades. Arrows show how typical points on curves shift 
for dual rotation of eight blades (four front and four back) 
Shifted points shown are for J = 5.1693. (See ref. 6.) 
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